"Barber pole turbulence" in large aspect ratio Taylor- Couette flow 
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Investigations of counter-rotating Taylor- Couette flow (TCF) in the narrow gap limit are con- 
ducted in a very large aspect ratio apparatus. The phase diagram is presented and compared to that 
obtained by Andereck et al. pi . The spiral turbulence regime is studied by varying both internal and 
external Reynolds numbers. Spiral turbulence is shown to emerge from the fully turbulent regime 
via a continuous transition appearing first as a modulated turbulent state, which eventually relaxes 
locally to the laminar flow. The connection with the intermittent regimes of the plane Couette flow 
(pCf) is discussed. 
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The "barber pole structure of turbulence" Q between 
two counter-rotating cylinders, also called spiral turbu- 
lence, is commonly described as alternating helical stripes 
of laminar and turbulent flow. There are few quantita- 
tive studies of this puzzling regime, where long range or- 
der coexists with small scale turbulence. In early studies 
Coles and Van Atta g-|^] measured the spiral shape in 
the mid-plane perpendicular to the cylinder axis and the 
pattern rotation rate. Later Andereck et al. reported 
qualitative observations when varying Ri [||. Hegseth 
et al. [|[0] described spiral turbulence within the frame- 
work of phase dynamics. All these studies were limited 
by their relatively small size. Only one helical turbulent 
stripe, winding no more than twice along the cylinder 
axis, could be observed. Altogether, the origin of this 
flow pattern remains unknown. 

Performing measurements in large aspect ratio Taylor- 
Couette flow, we show that the spiral turbulence bifur- 
cates continuously from the turbulent flow, appearing as 
a modulated turbulent state. After a rapid description of 
the experimental set up, we present the phase diagram 
and compare it to the one obtained by Andereck et al. |||] 
with a different cylinder radius ratio. Then we describe 
the successive steps leading to the fully turbulent flow 
before discussing the origin of spiral turbulence. Finally, 
we examine its breakdown into a spatio-temporal disor- 
dered regime similar to the laminar-turbulent coexistence 
dynamics observed in plane Couette flow HJ^. 

Let us recall here the characteristics of the experimen- 
tal set up |l^]. Our TC apparatus, two coaxial cylinders 
of radius n = 49.09 ± 0.005 mm and Vo = 49.95 ± 0.005 
mm with a useful length L = 380 ±0.1 mm, rotating inde- 
pendently, has a gap size d = 0.863 ±0.01 mm, large as- 
pect ratio Tz = L/d = 442 and Tq = 7r(ri ± ro)/d = 362, 
and a radius ratio 77 = 0.983 very close to 1. In this paper, 
length and time units are d and d'^ /u. r] being fixed, the 
flow is governed by the inner and outer Reynolds num- 
bers Ri^o = ^i,o^i,od/^i with Vti^o the angular velocities, 
and v the kinematic viscosity. The accuracy of Ri^o is 
better than 3%. 



We visualize the flow by a "fluorescent lighting "tech- 
nique [|lO| developed for this study. The water flow is 
seeded with Kalliroscope AQ 1000 (6 x 30 x 0.07/im 
platelets). The inner cylinder is covered by a fluores- 
cent film and the entire apparatus is UV-lighted. The 
fluorescent film re-emits a uniform visible lighting, trans- 
mitted through the fluid layer: the more turbulent the 
flow, the brighter it appears. As the gap is very thin, 
the Kalliroscope concentration is increased up to 25% by 
volume to enhance the contrast. A rheological study has 
shown that the fluid remains Newtonian, so that the only 
impact is a viscosity increase up to z^ = 1.13 10~^ m^/s at 
20° C. The flow is thermalized by water circulation inside 
the inner cylinder. At thermal equilibrium the tempera- 
ture is uniform in space up to 0.1°C and does not vary 
more than 0.1°C/hour. Images and spatio-temporal dia- 
grams (temporal recording of one line along the cylinder 
axis) are recorded by a CCD camera. Two plane mirrors 
reflect the two thirds of the flow hidden to the camera so 
that the whole cylindrical flow can be reconstructed. 



Coles Q has stressed the variety of flow structures ob- 
served according to the path followed in parameter space. 
The experimental procedure followed here is the most 
frequent one. The outer cylinder is slowly accelerated 
to the desired angular velocity while the inner one is 
at rest. Then, the inner cylinder is slowly accelerated 
{{1 / Q.i) .{d^i / dt) = 0.03% per second) in order to ensure 
quasi- static transitions. (It was checked that the Taylor 
vortices threshold is determined with an accuracy bet- 
ter than 1% as soon as the accelaration rate is less than 
0.05%). Recording the various flow regimes in parame- 
ter space (i?i, Ro)^ gives the experimental phase diagram 
(fig. |l|(a)) to be compared to that given by Andereck 
et al. [Q for a different geometry (77 = 0.878 and much 
smaller aspect ratio). The same regimes are observed but 
their locations in the parameter space change. The linear 
thresholds increase (at i?o = 0, the Taylor vortices appear 
for R^ = 320 instead of R* = 120). This was expected in 



regards of the plane Couette flow limit and agrees very 
wen with Snyder's law R* = 27[(1 - {r]))r]]^ - 5/3) [0. 
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(b) 
FIG. 1. Experimental phase diagram (77 = 0.983). (b) 
displays a zoom of the dotted rectangle of (a). The labels 
stand for AZI : azimuthal flow; TVF: Taylor vortex flow; 
WVF: wavy vortex flow; SPI&IPS: spiral and interpenetrat- 
ing spiral vortices; WIS: wavy interpenetrating spiral vortices; 
INT: intermittency; SPT: spiral turbulence; TUR: turbulence. 
The arrow indicates the Ro value below which the AZI flow 
jumps directly to INT (see text for details). The solid straight 
lines in (b), show the paths along which measurements are 
conducted (A: Ro = -922; B: Ro = -713; C: Ro = -641; 
D: Ri = 704). 

Accordingly, the other vortex flows (WVF, SPI, IPS...) 
also appear at higher Ri values. Second, the laminar- 
turbulent coexistence regions (INT and SPT), extend to 
higher Ro values. Finally, for Ro < —725, the azimuthal 
flow jumps directly to INT whereas, for Ro > —725, it 
bifurcates first to IPS or SPI. 

Figure g, shows snapshots of the whole flow, along 
paths A and C. For Ro = —922, the boundaries trig- 
ger turbulent stripes are at Ri lower than the expected 
linear instability threshold (fig. ||(a)). In fig. ||(b) peri- 
odically organized turbulent stripes have formed: spiral 
turbulence is established as a regular rotating pattern 
with uniform pitch angle along the cylinder axis. 




(a) 

FIG. 2. Turbulent spots and stripes along paths A and 
C. Left column is path A, Ro = -922: Ri = (a) 580, 
(b) 600, (c) 697, (d) 784. Right column is path C, 
R^ = -641: Ri = (e) 649, (f) 668, (g) 695, (h) 717. Each 
picture displays a 360 ° view of the whole flow (38 cm high 
and 31.4 cm wide. 



Further increasing i?^, the pattern maintains itself and 
the number of stripes increases (fig. ||(c)). Finahy, the 
pattern disappears smoothly as homogeneous turbulence 
takes over (fig. §|(d)). In this regime, one notices sev- 
eral domains of opposite pitch spirals. For Rq = —641, 
the azimuthal flow first bifurcates to a vortex fiow, within 
which, further increasing i?^, turbulent spots follow an in- 
termittent bursting dynamics (fig. ^e)). In fig. ^(f), the 
spots grow, merge and build up the turbulent spiral pat- 
tern. For illustration, we show, in fig. ||(g), two opposite 
pitch spirals competing and leading to a so-called "V"- 
shaped pattern. In figure 0(h), spiral turbulence turns 
into uniform turbulence as for Ro = —922. Fig. g dis- 
plays the spatio-temporal dynamics of a front between 
two opposite spirals. This front, moving along the cylin- 
der axis, is eventually absorbed at one of the extremities. 
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FIG. 3. Spatio-temporal diagram of spiral turbulence dis- 
playing a front between two domains of opposite pitch angle. 
Ro = —838, Ri — 707. The z-axis display the whole flow 
height (38 cm) and the total recording length is 20.48 s. 

Let us now consider the previous scenario in reverse. 
Decreasing Ri^ spiral turbulence appears from homoge- 
neous turbulence. At Ri = Ric^ a low amplitude mod- 
ulation, with wave vector k and frequency cj, uniformly 
emerges from the turbulent flow with no prefered direc- 
tion. A slow front dynamics occurs between domains of 
opposite pitch until the z ^ —z symmetry breaks and 
selects a single uniform spiral. Further decreasing Ri^ 
the modulation amplitude increases until a more com- 
plex dynamics occurs. The light intensity, which in the 
flrst approximation corresponds to the turbulence level, 
is recorded in spatio-temporal diagrams such as flg. H. In 
the following, we obtain the amplitude A, the axial com- 
ponent kz of the wave vector and the frequency co of the 
light intensity modulation by demodulation [g^. Owing 
to the azimuthal periodicity, the azimuthal component 
of the wavevector, kg = 2.ne/{ri + To), is quantized, ne 
being the number of pattern periods along one perimeter. 

Four paths in the {Ro^Ri) plan are investigated. We 
flrst consider the three vertical paths A, B and C {Ro = 
-922, Ro = -713 and Ro = -641). 
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FIG. 4. The spiral turbulence amplitude (a) and vertical 
wave number (b) vs. e. (Ro — —641 (*), Ro — —713 (o), 
R^ = -900 (•)). 

Fig. 11(a) displays the amplitude of the wave against the 
threshold distance e = \Ri — Ric\ / Ric- Ri is determined 
by extrapolation, assuming a linear dependence of A^ 
on Ri^ as suggested by the general framework of am- 
plitude equations. The resulting values, Ric = 820 for 
Ro = -922, Ric = 785 for Ro = -713 and Ric = 740 
for Ro = —641, agree with the observed thresholds. The 
axial component kz of the wavevector decreases linearly 
with e from a critical value /c^ close to 0.30 (see flg. §(b)). 
As long as the spiral is well deflned, the azimuthal wave 
number no equals 6 along paths A and B and 5 along 
path C. no being constant and kz decreasing with e, the 
spiral pitch, (p = diTctdin{ko/kz)^ increases with e from a 
critical value 0"^ = 18° up to 40°. 
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FIG. 5. The angular velocity cOs vs. Qm, (Ro = —641 (*), 
R^ = -713 (o), Ro = -900 (•) and Ri = 704 (+)). 

We now consider the horizontal path D with Ri = 704. 
For Ro > —630, the flow is fully turbulent. Decreasing 
Ro^ spiral turbulence emerges from the fully turbulent 
regime. The initial azimuthal wave number no = ^ in- 
creases to 7161 = 6 by Eckhaus-like re-arrangements. The 
axial wave vector component remains constant on the ex- 
plored values —950 < Ro < —630. Finally, flg. ^displays 
COs = ^ I^Qi the angular velocity of the spiral turbulence 
pattern, vs. (^^ = \(Q.i^Vlo)- As seen by Coles and Van 
Atta ||-||] , the pattern rotates with the mean angular ve- 
locity of the two cylinders (cj^ = Q^rn) and not with the 
outside cylinder as reported by Andereck et al. [|l|] as well 



as by Goharadzeh and Mutabazi [Q . In the latter study, 
the range of variation for Vti/Vto is actually too small to 
distinguish a dependence on VLq from one on Q^rn- 

Altogether, the spiral turbulence emerges as a wavy 
turbulent state pattern, which bifurcates continuously 
from the homogeneous turbulent flow. This wavy be- 
havior agrees with the idea of a phase dynamics analy- 
sis suggested by Hegseth et al. [||J^] to describe the pitch 
variation, they observed. However the present study does 
not report any pitch variation. In their analysis, Hegseth 
et al, on the basis of the mean flow measurements by 
Coles and Van Atta p], suggest that the spiral turbu- 
lence induces a backflow, whose mean axial component 
is differently modified at both extremities. The resulting 
different pitches locally imposed in the boundary con- 
ditions drive the pitch variation along the axis. In the 
present study, both the large aspect ratio and the gap 
thinness considerably reduce the backflow and the influ- 
ence of the end boundaries. 

For low values of Ri , we observed disconnected turbu- 
lent domains surrounded by a laminar flow, which evokes 
the transitional dynamics observed in the plane Couette 
flow [Q, where turbulent domains move, grow, decay, 
split and merge leading to the so-called spatio-temporal 
intermittency (STI), a process in which active/turbulent 
regions may invade absorbing/laminar domains out of 
which disorder cannot spontaneously emerged |1^ . Con- 
sidering the vicinity of the plane Couette flow limit, 
?7 = 1, and following Colovas and Andereck llg], one 
would like to know whether Taylor- Couette flow actually 
exhibits STI and how the flow evolves from a continuous 
turbulent state modulated in space and time to two-state 
coexistence. In a first attempt to answer these questions, 
we binarized the flow into alternating laminar and tur- 
bulent domains as is common in an STI context. This 
first-order approximation should not influence the results 
presented below, at least at a qualitative level. 
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FIG. 6. A (Ro = -713 (o) and Ro = -922 (•)), the pattern 
wavelength and C (Ro = -713 (A) and Ro = -922 (*)), the 
turbulent stripes width, against e. 

Fig. 1^ displays the turbulent stripes width, (^ and the 
wavelength, A of the spiral pattern, parallel to the wave 
vector direction. A increases faster than (: the minima of 



the modulation becomes more and more distant from the 
most active turbulent regions. Simultaneously the modu- 
lation amplitude increases and these minima are less and 
less turbulent. As a result they become areas favorable 
to local relaxation to the laminar flow and the modulated 
turbulent state turns into a periodic pattern of alternat- 
ing turbulent and laminar stripes separated by fronts. 
This totally different spatio-temporal dynamics rapidly 
turns into the disordered STI regime. 

Our very large aspect ratio TC apparatus has allowed, 
for the first time, to visualize and quantitatively analyse 
the "barber pole turbulence" at large scales. First, this 
study brings quantitative evidence of the spatio-temporal 
periodicity of the pattern and its dependance on control 
parameters. Second, it shows that the spiral turbulence, 
previously described as alternating turbulent and lam- 
inar bands, actually appears through a modulation of 
the homogeneous turbulent state. To our knowledge, it 
is the first time that the transition from homogeneous 
turbulence to spiral turbulence is investigated. Finally 
measurements of the turbulent stripes width relatively to 
the pattern wavelength have allowed us to describe the 
relaxation mechanism followed between the modulated 
turbulent state and the laminar-turbulent coexistence. 
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technical assistance. 
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